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ABSTRACT  OF  OBJECTIVES  AND  ACCOMPLISHMENTS 

The  principal  focus  of  this  research  program  is  the  electronic  structure  aspects  of 
nonadiabatic  processes  related  to  the  stability,  formation  and  detection  of  high  energy  density 
materials(HEDMs).  During  the  current  grant  our  work  has  focused  on  metal  doped  hydrogen  as 
cyrogenic  fuels.  These  studies  were  enabled  by  a  unique  system  of  computer  codes,  the 
BROOKLYN  programs,  that  provide  advanced  capabilities  for  the  study  of  the  electronic  structure 
aspects  of  spin-conserving  and  spin-nonconserving  electronically  nonadiabatic  processes  and  have 
permitted  us  to  make  significant  contributions  in  this  area  of  chemistry.  Our  most  significant 
algorithmic  accomplishment  in  this  grant  period  has  been  the  development  of  algorithms  for 
locating  intersections  of  distinct  seams  of  conical  intersections.  These  algorithms  proved 


invaluable  in  analyzing  the  seams  of  conical  intersection  in  AI+H2,  B+H2  and  C+H2  which  are 
relevant  to  the  kinetic  stability  of  the  corresponding  van  der  Waals  species  in  cryogenic  H2.  We 
have  also  investigated  the  stability,  and  lifetime  with  respect  to  spin-nonconserving  radiationless 
decay  of  the  C2A  state  of  BAr  and  BNe  and  the  B3n  state  of  CAr.  These  studies  of  metal  rare  gas 
bonding  support  the  experimental  work  of  P.  J.  Dagdigian  on  the  detection  of  metal  -rare  gas 
complexes  and  the  work  of  P.  Langhoff  on  the  spectral  theory. 
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REPORT  DEFINITION 


F49620-96- 1-0017 
David.  R.  Yarkony 

On  the  Origin  of  Metastability  in  Energetic 
Species 

Final  Technical  Report 
(15  December  1995  -14  December  1998) 

SUMMARY 

This  grant  period  was  one  of  considerable  excitement  and  discovery.  Our  studies  of  A1  + 
H2  led  us  to  explore  novel  topologies  for  seams  of  conical  intersections  and  ultimately  to  develop 
an  algorithm  for  identifying  intersecting  seams  of  conical  intersections,  our  most  significant 
algorithmic  accomplishment  in  this  grant  period.  This  algorithm  proved  invaluable  in  analyzing  the 
seams  of  conical  intersection  in  B+H2,  C+H2  and  Li+  H2  which  are  relevant  to  the  kinetic  stability 
of  the  corresponding  van  der  Waals  species  in  cryogenic  H2.  Our  C+H2  study  is  particularly 
noteworthy  since  it  clearly  illustrated  the  potential  pitfalls  in  characterizing  seams  of  conical 
intersection  without  the  above  noted  algorithm.  We  have  also  investigated  the  stability  the  12A 
state  of  B Ar  and  of  BNe  and  the  B3n  state  of  CAr,  and  determined  their  lifetimes  with  respect  to 
spin-nonconserving  radiationless  decay.  These  studies  of  metal-rare  gas  bonding  and  radiationless 
decay  support  the  experimental  work  of  P.  J.  Dagdigian  on  the  detection  of  metal  -rare  gas 
complexes  and  the  work  of  P.  Langhoff  on  the  spectral  theory. 
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TECHNICAL  REPORT 


The  thrust  of  our  research  program  is  two  fold:  (i)  to  develop  computational  techniques  to 
determine  the  stability  with  respect  to  radiationless  decay  of  potential  high  energy  density  materials 
and  (ii)  to  use  of  those  techniques  to  study  chemically  relevant  systems  and  assess  the  feasibility 
of  their  spectroscopic  characterzation/detection.  During  the  current  grant  have  developed  new  tools 
for  the  study  of  nonadiabatic  processes.  In  addition  to  their  use  in  accessing  the  stability  of 
energetic  materials  these  algorithms  have  potentially  broad  impact  in  the  field  of  nonadiabatic 
chemistry.  Below  we  discuss  the  most  significant  components  of  our  current  research  effort. 

ACCOMPLISHMENTS/NEW  FINDINGS 

I.  Characterizing  novel/unexpected  seams  of  conical  intersection 

As  part  of  our  Air  Force  sponsored  research  we  have  used  the  criterion1  gIJ  (Rx)  x 
h^Rx)  =tIJ  (Rx)=  0  where  Rx  is  a  point  on  a  C2v  seam  of  conical  intersection  to  establish  the 
existence  of  a  Cs  seam  intersecting  the  C2v  seam  at  Rx.  Here 


2g"(R) 

S  (c'(Rf )  - c'(R,))f ^pV(R ,)  +  c' (Rxj) 
s  c'(R,)^(R)c'(R,X 

(la) 

fc"(R) 

(lb) 

E/(R)]c/(R)  =  0 

(lc) 

H{ R)  is  the  electronic  hamiltonian  matrix,  and  Ej(Rx)  =Ej(Rx).  It  is  important  to  note  that  g^(Rx) 
x  h^Rx)  is  readily  evaluated  using  analytic  gradient  techniques  in  algorithms  we  have  developed. 
This  approach  will  be  particularly  valuable  in  general  polyatomic  molecules  ( >  3  atoms )  where  the 
seam  of  conical  intersection  is  a  surface.  In  such  situations,  in  the  absence  of  this  criterion,  it 
would  be  difficult  to  distinquish  between  intersecting  seams  and  the  dimensionally  allowed 
continuation  of  a  the  seam  into  regions  of  reduced  symmetry. 

H:  Conical  Intersections  and  Energetics  Materials 
(A)  C(3P)  +  H2,  B(2P)+  H2,  A1(2P)+  H2 

The  reactions 

M(2S+lp)+H2-4  MH2(l2S+lAi,  l2S+lBi)  — »  MH2(2S+1  Y)  (2a) 

MH(X!Z+)  +  H  (2b) 
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where  (2S+1,  i,  i',Y)  =  (2,l,2,Ai)  for  aluminum3  and  boron4  and  (2S+1,  i,  i',  Y)  =  (3,2,1,  Bj) 
for  carbon,5  are  relevant  to  the  stability  and  properties  of  M  doped  cryogenic  hydrogen.  The 
formation  of  the  dihydride  MH2(2S+1Y),  channel  (2a),  could  limit  the  stability  of  the  van  der 
Waals  complex  M  —  H2  that  constitutes  the  energetic  material,  while  channel  (2b)  may  be 
involved  in  the  combustion  of  the  energetic  material. 

Low  energy  pathways  to  the  dihydride(a  2Ai  state  for  aluminum  and  boron  and  a  3B2  state 
for  carbon)  involve  the  2Bi  section  of  the  12A'  PES  for  aluminum  and  boron  and  the  3A2  section 
of  the  13A"  PES  for  carbon.  Thus  the  C2V  constrained  minimum  energy  pathway  to  the 
dihydride,  T-shaped  approach  of  M  to  H2,  precedes  over  the  ridge  corresponding  to  the  symmetry- 
allowed  2S+1Bi  -  2S+1Ai-  seam  of  conical  intersection.  The  minimum  energy  crossing  point 
(MECP),  the  minimum  energy  point  on  the  2S+lBi  -  2S+1  A,'  seam,  is  the  transition  state  for  this 
path.  However  there  are  two  directions  perpendicular  to  the  ridge  one  that  preserves  C2v 
symmetry,  the  aforementioned  constrained  reaction  path,  and  an  asymmetric  mode  that  removes  the 
C2v  symmetry.  Displacements  along  this  mode  give  rise  to  general  Cs  geometries  which 
conventionally,  would  be  devoid  of  conical  intersections,  and  would  include  the  true  transition 
state.  In  fact  the  transition  state  for  these  systems  is  found  quite  near  the  MECP  which  therefore 
provides  a  useful  upper  bound  to  the  barrier  for  forming  the  dihydride. 

The  MECP  for  the  BH2  12A'  -  22A'  conical  intersection  seam  was  found  to  be 
approximately  0.6eVabove  the  B  +  H2  asymptote.  This  result,  which  is  less  than  half  the  value  in 
AIH2,  represents  a  considerable  lowering  of  a  previous  estimate6  of  the  insertion  barrier  and  thus 
has  potentially  important  implications  for  the  stability  of  B  doped  cryogenic  hydrogen.  For  CH2 
the  MECP  on  the  3A2  -  3B 1  C2V  seam  of  conical  intersection  was  found  to  be  only  0.77  kcal/mol 
above  the  C(3P)  +  H2  asymptote.  Thus  the  barrier  for  the  van  der  Waals  complex  — »  dihydride 
conversion  is  at  best  small. 

These  results  are  of  considerable  practical  importance.  Yet  they  may  not  be  the  most 
influential  results  of  these  study.  The  figure  below7  presents  a  schematic  representation  of  the 
1,22A’  seam  in  A1(2P)H2  reporting  Rx(r)  =  (R(r),  r,  y (r))  for  which  El2A.(Rx(/'))  =  E22A,(Rx(r)), 
that  is  the  locus  of  points  on  the  12A’  -  22A'  seam  of  conical  intersection  using  Jacobi  coodinates, 
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R  =  ( R ,  r,  f)  with  r  as  the  seam  parameter.  It  would  be  expected  that  yO)  =  90°  for  all  r,  that  is 
the  seam  should  have  C2V  symmetry.  However  this  is  plainly  not  the  case!  For  r  <  ~  3.  lao,  y(r) 
=  90°  as  expected. 


However  for  each  r  >  ~3.1ao  there  are  three  points  of  conical  intersection.  One  point  has  y  =  90° 
and  is  part  of  the  C2v  seam.  There  are  also  two  equivalent  points  with  /  =  yx  ,  180°  —  yx  having 
only  Cs  symmetry.  It  is  important  to  emphasize  that  all  R*  are  on  the  12A'  -  22A'  seam  of  conical 
intersection.  This  unusual  feature  represents  a  trifurcation  of  the  C2V  seam  as  r  increases  past 
~3.1ao.  Equivalently  it  is  the  intersection  of  a  Cs  seam  and  a  C2v  seam  at  the  C2v  point  Rdd-  Since 
points  of  conical  intersection  are  referred  to  as  diabolical  points,^  points  at  the  intersection  of  two 
seams  of  conical  intersection  are  referred  to  as  doubly  diabolical  points.1 

The  importance  of  this  feature  cannot  be  overemphasized.  Its  existence  means  that  the 
significant  nonadiabatic  effects  associated  with  conical  intersections  cannot  automatically  be 
assumed  to  exist  only  for  nuclear  configurations  yielding  symmetry-allowed  intersections. 
Additional  confluences  are  possible  in  unexpected  regions  of  nuclear  coordinate  space  and  may 
result  in  'confluences  of  confluences’. 
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In  view  of  the  A1(2P)-H2  result  it  was  natural  to  ask  whether  a  similar  result  would  be 
obtained  in  the  B-H2  and  C  -  H2  systems.  The  Cs  seam  in  AIH2  was  located  quite  serendipitously. 
The  criterion  described  in  section  I  allowed  us  to  anticipate  the  existence  of  the  Cs  seam  given  only 
knowledge  of  the  (standard)  C2V  seam.  The  following  two  figures4  consider  the  BH2  12A'-22A' 
seam  of  conical  intersection.  The  first  figure  presents  Rx(r)  =  (R(r),  r,  y(r)  =  90)  and  £l2A,(Rx(r)) 
=  E22  a  ,(Rx(r))  relative  to  B(2P)  +  H2.  This  is  a  standard  C2V  seam  of  conical  intersection  between 
a  2Ai  and  a  2B2  state. 


However  as  the  following  figure  shows  gIJ  x  vanishes  near  r  =  3.1ao  (*  point  above) 
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1.0 . II  ■  I  .  T  .  I 
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Thus  an  intersecting  Cs  branch  of  the  seam  is  expected  at  R*(r~3.1).  The  following  figure 
presents  that  intersecting  branch,  Rx($  =  (R(y),  KY)>  Y)  and  ^a'^xW)  =  ^A'^xC)^)- 


Y(deg) 
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The  above  result,  published  as  a  communication  in  the  Journal  of  Chemical  Physics,  demonstrates 
the  power  of  this  approach  as  a  general  tool  for  locating  intersecting  seams  of  conical  intersection.4 

While  these  results  were  quite  unexpected  perhaps  the  most  fundamentally  surprising  result 
to  emerge  these  studies  was  the  demonstration  that  for  CH2  in  addition  to  the  conventional  C2V 
and  CooV  seams  of  conical  intersection,  there  again  exists  a  Cs  seam  of  conical  intersection  that  in 
turn  intersects  the  C2V  seam.  This  result  was  unexpected  since  for  CH2  the  C2V  (  and  CooV )  seams 
of  conical  intersection  had  been  carefully  studied  by  Harding,  Schatz  and  coworkers9’10  as  part  of 
a  determination  of  global  1,23A"  PESs  based  on  6000  high  quality  ab  initio  points!  For  this 
reason  the  CH2  results  appeared  as  a  communication  in  J.  Chem.  Phys.5  The  CH2  results  are 
expected  to  be  particularly  important  since  the  MECP  on  the  Cs  portion  of  the  seam  as  exoergic 
relative  to  the  CH(X2n)  +  H(2S)  asymptote. 

The  trifurcation  exhibited  by  CH2,  BH2  and  AIH2  is  highly  nonstandard  it  was  important 
to  ask  how  it  could  be  demonstrated  that  the  Cs  seam  of  conical  intersection  represents  a  true 
intersection  seam  rather  than  narrowly  avoided  crossings.  One  could  determine  whether  in  fact 
these  points  are  points  of  conical  intersection  by  considering  the  phase  of  the  electronic 
wavefunctions  for  closed  paths  C  surrounding  the  point.  However  this  can  be  quite  tedious.  We 
have  shown11  that  X(Cp),  the  circulation  of  the  derivative  coupling,  fIJ(R),  where 
frJ,( R)  =  ('F/(r;R)|(<9/<?T)'F/(r;R))r  and  x  is  an  internal  coordinate,  along  a  small  loop  of  radius  p 
around  the  point  in  question,  Cp,  can  be  used  to  consider  this  question  since: 

X(Cp)=ffw(R)-^R-^x7r  (3) 

cp 

where  K  =  0(1)  if  0(1)  conical  intersections  are  enclosed  inCp.  This  approach  has  the  advantage 
that  the  phase  of  the  integrand  at  the  endpoint  of  the  loop  is  known,  since  f^(R)  does  not  change 
sign  after  traversing  a  closed  loop.  This  provides  a  useful  control  in  assigning  the  phase  of  the 
integrand  at  neighboring  points.  The  existence  of  three  independent  components  for  the  derivative 
coupling  vector  fIJ(R)  can  also  help  to  decide  phase  relationships  in  otherwise  ambiguous 
situations  without  the  need  to  determine  additional  points.  The  use  of  eq.  (3)  is  made 
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computationally  tractable  by  an  algorithm  we  have  developed12,13  for  the  efficient  evaluation  of 
the  fIJ  using  analytic  gradient  techniques. 

Eq.  (3)  was  used  to  confirm  that  for  AIH2,  BH2  and  CH2  the  1,22S+1A  potential  energy 
surfaces  do  in  fact  exhibit  this  remarkable,  and  heretofore  virtually  unknown,  locus  of  seams  of 
conical  intersection  (the  only  other  reports  based  on  ab  initio  wavefunctions  are  for  the  1]A'  and 
21  A'  states  of  O314  and  the  23A"  and  33A"  states  of  CH215). 

(B)Li  +  H216 

The  ultimate  fate  of  excited  Li(2P)  in  cryogenic  H2  matrices  is  relevant  to  the  quantitative 
detection  of  ground  state  Li(2S)  in  this  matrix.  The  quenching  of  electronically  excited  alkali  metal 
atoms  has  been  found  by  Scoles  and  coworkers17  to  be  extremely  efficient.  Fajardo's  group  has 
also  studied  this  quenching  process.  We  have  investigated  the  12A'  -  22A’  conical  intersection 
seam  in  UH2.  These  states  correlate  with  Li(2S)+H2  and  Li(2P)+H2  so  that  this  seam  is 
principally  responsible  for  this  quenching.18,16  We  have  again  demonstrated  the  existence  of  a 
previously  unknown  portion  of  this  seam  with  Cs  symmetry  that  intersects  the  C2V  seam.  This 
results  is  quite  surprising  in  view  of  the  long  history  of  theoretical  study  of  these  states  (see 
references  in  Ref.  16). 

HI:  Spin-nonconserving  Radiationless  Decay  and  its  relation  to  Laser  Detection  of  HEDMs 
A.Unusual  Bonding  and  the  heavy  atom  spin-orbit  effect  in  the  open  shell  van  der  Waals  BAr 

The  importance  of  boron  doped  in  cryogenic  hydrogen  matrices  to  the  cryogenic  fuels 
program  lead  Dagdigian  and  coworkers  to  systematically  study  the  interactions  of  the  open  shell 
van  der  Waals  molecules  formed  by  boron  in  its  electronic  ground,  or  excited,  state  with  a  rare  gas 
or  molecular  hydrogen.1^'23  That  research  suggested  that  the  C2A  state  of  BAr  is  quite  strongly 
bound  but  is  rapidly  predissociated  by  spin-orbit  induced  coupling  to  the  repulsive  14II  state. 
The  later  observation  seemed  counterintuitive  since  spectroscopic  studies  suggested  the  open  shell 
moiety  to  be  a  largely  unperturbed  boron  atom  which  has  at  best  modest  spin-orbit  interactions. 
For  this  reason  we  considered  the  electronic  structure  of  the  C2A  and  l4n  states  of  BAr  and  the 
spin-orbit  interactions  that  couple  these  states.  The  results  of  this  study24  were  quite  surprising 
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or 

and  had  important  implications  for  a  long  standing  controversy  concerning  the  anomalous  fine 
structure  splitting  in  alkali  metal  -  rare  gas  dimers26  and  for  other  AFOSR  research  efforts, 

'll 

including  M.  Duncan's  work  at  the  University  of  Georgia  on  metal  ion-rare  gas  interactions  and 
P.  W.  Langhoff  s  work  on  spectral  theory  28 

The  computational  treatment  was  particularly  challenging  owing  to  the  need  to  treat  near¬ 
degeneracy  effects  attributable  to  the  Rydberg  ls22s23d]  2D  electron  configuration.  The  results  of 
the  electronic  structure  treatment,  which  are  pictured  below,  confirm[C2A(ab)]24  the  unusually 
strong  binding  in  the  BAr  C2A  electronic  state  observed  by  Dagdigian[C2A(exp)]. 
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More  significantly  however  it  emerged  that  the  predissociation  of  the  C2A  state29  was  a 
consequence  of  the  heavy  atom  effect  on  the  C2A  -  14IT  spin-orbit  coupling,  HS0(C2A,  l4n).  The 
large  increase  in  Hso(C2A,  14IT)  near  the  C2A  equilibrium  geometry  seen  in  the  above  figure 
represents  the  heavy  atom  effect,  the  borrowing  of  spin-orbit  coupling  by  boron  from  argon.  The 
marked  increase  of  Hso(C2A,  l4n)  with  decreasing  internuclear  separation  suggested 
nonnegligible  overlap  between  valence  orbitals  on  argon  and  boron  which  lead  us  to  establish  a 

novel  origin  for  the  bonding  in  the  C2A  state  -  see  below. 

(B)  Extension  to  BNe(C2n) 

We  extended  this  work  to  the  C2A  state  of  BNe.  Dagdigian  and  coworkersrecorded  the 
BAr23  and  BNe3^  laser  fluorescence  excitation  spectra  and  the  BAr  fluorescence  depletion  spectra 
in  the  region  of  the  atomic  boron  2s2p2  2D  <—  2s2p  2P.  The  C2A  state  of  BNe  was  found  to 
fluoresce  but  that  of  BAr  did  not.  The  absence  of  radiative  decay  in  BNe(C2A)  suggests 
predissociation  by  a  repulsive  4n  state.  To  help  interpret  these  results  we  performed  computational 
studies  to  determine  the  mechanism,  and  rates,  of  BNe(C2A,  v)  predissociation.  The  potential 
energy  curves  for  both  the  C^A  state  and  the  repulsive  l4n  state,  and  the  spin-orbit  interactions 
that  couple  these  states  were  determined  and  used  to  compute  the  radiationless  decay  rates. 
Combined  with  the  above  BAr  results  these  calculations  provide  a  conceptually  clear  picture,  and 
quantitively  reliable  treatment,  of  the  BRg  (C2A,v)  predissociation.29’3 

(C)  Implications 

While  the  BRg  (C2A,v)  predissociation  rates  are  of  practical  significance  to  the  the  HEDM 
program  perhaps  the  most  important  outcome  of  these  studies  are  the  generalizable  insights 
concerning  the  connection  between  the  bonding  in  these  species  and  the  heavy  atom  effect,  the  later 
being  principally  responsible  for  the  BRg  (C2A,v)  predissociation.  The  orbital  overlap  that  leads  to 
the  heavy-atom  effect  in  the  spin-orbit  coupling  l4n  ~  C-A.  gives  rise  to  overlap  of  a  vacant  axial 
M  orbital  and  the  highest  occupied  axial  orbital  on  Rg.  Pursuing  this  line  of  reasoning  we  showed 
that  the  strong  binding  is  best  described  as  a  dative,  or  coordinate  covalent,  bond  with  the  Ar 
furnishing  the  electrons  from  its  fully  occupied  3pz  orbital  to  an  empty  B  2s2pz  hybrid  orbital. 
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Overall  the  binding  results  from  a  balance  of  this  bonding  interaction  and  electron  correlation 
effects,  with  high-order  configuration  interaction  effects  being  very  important  in  describing  the 
overall  binding.  Since  a  spectral  theory  treatment  of  bonding  in  metal-rare  gas  van  der  Waals 
complexes  depends  on  accurate  knowledge  of  the  pairwise  interactions  of  both  the  ground  and 
excited  electronic  states28  this  finding  has  important  implications  for  P.  W.  Langhoffs  spectral 
representation  work. 

On  the  basis  of  this  analysis  we  were  able  to  explain26  the  fine  structure  splitting  in  the 
2F1q  states  of  Li(2p)Ne  and  Li(2p)Ar  ,  which  had  been  a  matter  of  considerable  debate.25  It  had 
been  previously  been  asserted25  that  mixing  of  Rydberg  states  on  the  rare  gas  with  metal  valence 
states  accounted  for  the  observed  heavy  atom  effect.  However  our  model  of  the  bonding  indicated 
that  valence  orbitals  on  the  rare  gas  should  be  responsible  for  the  observed  effect.  We  were 
subsequently  able  to  show26  that  the  heavy  atom  effect  was  in  fact  attributable  to  valence  orbitals 
on  the  rare  gas.  These  notions  are  also  relevant  to  the  fine  structure  splittings  in  eg.  the  MgRg  and 
MgN2  complexes  studied  by  Professor  M.  Duncan  as  part  of  his  AFOSR  sponsored  work  at  the 
University  of  Georgia. 

(iv)  CAr 

Finally  we  considered  the  B3n  state  of  CAr.32  This  state  is  potentially  of  use  in  a  laser 
induced  fluorescence  detection  scheme  for  the  van  der  Waals  complex  CAr(X3Ir).  It  is  however 
predissociated  by  the  5E-  state  correlating  with  C(ls22s2p3).  The  figure  below  summarizes  our 
electronic  structure  results. 
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The  B3n  state  is  quite  strongly  bound  for  a  van  der  Waals  complex,  with  binding  energy 
-5000cm-1 !  The  bonding  and  existence  of  a  clear  heavy  atom  effect  in  the  5Ir  -  B3n  spin-orbit 
interaction  are  consistent  with  the  dative  bonding  model  noted  above.  We  find  that  spin 
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nonconserving  predissociation  of  the  the  B3n  state  will  not  preclude  use  of  the  B3n  state  to  detect 
the  ground  X3X~  state  although  Franck-Condon  factors  limit  the  vibrational  levels  of  the  B3n  state 
that  are  accessible  in  practice. 

Future  Directions 

As  part  of  our  renewal  proposal  we  will  work  with  David  Weeks  of  AFTT  on  the 

B+H2  ->  BH  +H 

scattering  problem.  We  are  in  the  process  of  generating  the  1,22A'  potential  energy  surfaces  and 
the  derivatve  couplings.  An  approximate  transformation  to  diabatic  states  is  being  developed  'on 
the  fly'  based  on  our  analysis  which  shows  the  transformation  in  question  necessarily  removes  all 
the  singularity  in  the  derivative  coupling.33  It  is  important  to  observe  that  this  approach  is  capable 
of  handling  the  trifurcation  of  the  C2v  seam  of  conical  intersection  described  above.  Given  the 
novel  seam  of  conical  intersection  in  this  system  these  potential  enegy  surfaces  and  the  subsequent 
dynamical  studies  should  be  of  general  interest  for  the  field  of  nondiabatic  dynamics  as  well  as  of 
practical  concern  in  the  HEDM  program. 
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